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tion for the macroscopic gravitational potential is derived. This is a biharmonic 
equation which is a non-trivial generalization of the Poisson equation of Newtonian 
gravity. In case of the strong gravitational quadrupole polarization it essentially 
holds inside a macroscopic matter source. Outside the source the gravitational 
potential fades away exponentially. The equation is equivalent to a system of the 
Poisson equation and the nonhomogeneous modified Helmholtz equations. The gen- 
eral solution to this system is obtained by using Green's function method and it 
does not have a limit to Newtonian gravity. In case of the insignificant gravita- 
tional quadrupole polarization the equation for macroscopic gravitational potential 
becomes the Poisson equation with the matter density renormalized by the factor 
including the value of the quadrupole gravitational polarization of the source. The 
general solution to this equation obtained by using Green's function method has a 
limit to Newtonian gravity. 
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1 Introduction 



General relativity as a classical theory of gravity is known to have some remarkable analo- 
gies with classical Maxwell's macroscopic theory of electromagnetism (see, for example, pQ 
and references therein). The physical motivation and intuition for many of the problems 
posed in general relativity originate therefore in their electromagnetic analogies where 
physics and formalism are much easier to deal with. One of such problems of the pri- 
mary importance is that of gravitational waves which has been inspired and put forward 
mostly owing to our deep understanding of the structure and physics of electromagnetic 
waves. Some important issues in the physics of gravitational waves in general relativity, 
however, remain obscure. The questions whether or not gravitational waves undergo the 
refraction in a gravitating macroscopic (continuous) media, the speed of a gravitational 
wave changes (i.e. slows down) in a material medium, the phenomenon of gravitational 
polarization exists, have been not approached and even properly posed as yet. 

To determine the structure of the macroscopic energy-momentum tensor resulting from 
averaging out a microscopic matter source, the problem of construction of a continuous 
(macroscopic) matter model for a given point-like (microscopic) matter distribution in 
general relativity has been formulated in [2]. The existing approaches have considered 
and a physical analogy with the similar problem in classical macroscopic electrodynamics 
has been pointed out. The procedure due to Szekeres |3] in the linearized general relativ- 
ity on Minkowski background space-time to construct a tensor of gravitational quadruple 
polarization by applying Kaufman's method of molecular moments jl] for derivation of 
the polarization tensor in macroscopic electrodynamics and to derive an averaged field 
operator by utilizing an analogy between the linearized Bianchi identities and Maxwell 
equations, has been analyzed. The approach of Szekeres to construct a tensor of the grav- 
itational quadruple polarization is based on the following assumptions: (a) the linearized 
theory of gravity on Minkowski background space-time; (b) the linearized field equations 
are taken as the linearized Bianchi identity to employ an analogy between gravitation and 
electromagnetism; (c) the covariant method of molecular moments of Kaufman is applied 
to construct a tensor of quadruple gravitational polarization. The procedure is shown to 
possess some inconsistencies, in particular, (1) it has only provided the terms linear in 
perturbations for the averaged field operator which do not contribute into the dynamics 
of the averaged field, and (2) the analogy between electromagnetism and gravitation does 
break upon averaging. A macroscopic gravity approach in the perturbation theory up 
to the second order on a particular background space-time taken to be a smooth weak 
gravitational field has been applied to write down a system of macroscopic field equations 
|2] 0, M Isaacson's equations jTj with a source incorporating the quadruple gravita- 
tional polarization tensor, Isaacson's energy-momentum tensor of gravitational waves [S] 
and energy-momentum tensor of gravitational molecules and corresponding equations of 
motion. The system of equations is shown to be underdetermined. A suitable set of 
material relations which relate all the tensors has been proposed [2|, [5], jB], so that the 
full system of the field equations and the material relations become determined. 

In this paper the system of equations is used to find a solution to the Szekeres model of 
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the gravitational quadrupole polarization. A model for the static weak-field macroscopic 
medium is analyzed and the equation for the macroscopic gravitational potential is de- 
rived. This is a biharmonic equation which generalizes the Poisson equation of Newtonian 
gravity. In case of the strong gravitational quadrupole polarization it essentially holds 
inside a macroscopic matter source. Outside the source the gravitational potential fades 
away exponentially. The equation is equivalent to a system of the Poisson equation and 
the nonhomogeneous modified Helmholtz equations. The general solution to this system 
is obtained by using Green's function method. This solution does not have a limit to 
Newtonian gravity. In case of the insignificant gravitational quadrupole polarization the 
equation for macroscopic gravitational potential becomes the Poisson equation with the 
matter density renormalized by the factor including the value of the quadrupole gravita- 
tional polarization of the source. The general solution to this equation has been obtained 
by using Green's function method and it has a limit to Newtonian gravity. 

2 Macroscopic Description of Gravity 

To consider the problem of macroscopic description of gravitation, an approach of macro- 
scopic gravity has been proposed earlier in [5] , [TUJ (see [TTj , (T21 , P3I for discussion of the 
problem and references therein, jHj for discussion of the physical status of general rela- 
tivity as either a microscopic or macroscopic theory of gravity). A covariant space-time 
volume averaging procedure for tensor fields , [ID] , [13 has been defined and proved to 
exist on arbitrary Riemannian space-times with well-defined properties of the averages. 
Upon utilizing the averaging scheme, the macroscopic gravity approach has shown that (i) 
averaging out Cartan's structure equations brings about the structure equations for the 
averaged (macroscopic) non-Riemannian geometry and the definition and the properties 
of the correlation tensor, (m) the averaged Einstein's equations become then the macro- 
scopic field equations and they must be supplemented by a set of differential equations for 
the correlation tensor, (m) it is always possible to extract the field operator of the same 
form as that of Einstein's equations for the Riemannian macroscopic metric tensor and its 
Ricci tensor while all other non-Riemannian correlation terms go to the right-hand side 
of the averaged Einstein's equations to give geometric correction to the averaged (macro- 
scopic) energy-momentum tensor. It is been also shown [TI], [T2| , ^] that only in case 
of neglecting all correlations of the gravitational field the averaged equations becomes 
the macroscopic Einstein equations for the Riemannian macroscopic metric tensor and 
its Ricci tensor with a continuous matter distribution. This result reveals the physical 
status of using the so-called standard procedure in cosmology when one claims that the 
Einstein equations preserve their form after substitution of a discrete matter model by a 
continuous one. 

In using the Einstein equations for a matter distribution in the form of a set of point- 
like mass constituents, there is a problem of adequate application, or validity, of the 
Einstein equations when such a matter distribution is substituted by a continuous matter 
distribution while the field operator in the left-hand side of the equations is kept un- 
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changed. This problem as it stands in cosmology 1 is called the averaging problem Pj, 
[T7j . [TBI, [13, EH!- Indeed, let us consider the Einstein equations in the mixed form 2 

9 ae r t , - i($rv = -^ (discrete) (1) 

with 

tf^^^t^x-z^)} (2) 

i 

where t^" is a energy-momentum tensor for a point-like mass moving along its world line 
z^ = z^(T(i)) parameterized by T(j) and % counts for the matter particles in the distribution 
(J2J). Changing the discrete matter distribution to a continuous (hydrodynamic) one in the 
right-hand side of (JTJ), which is the standard approach in cosmology 0, [Hj, [H], [H], EH] 
made phenomenologically on the basis of assumption about the uniformity and isotropicity 
of distribution of galaxies, or cluster of galaxies, throughout the whole Universe, means 
an implicit averaging denoted here by (•) 

^(discrete) (x) _ {%) = (^t^X - 3 (0 (t (1) )]\ . (3) 

Given a covariant averaging procedure (•) for tensors on space-time, the averaging out of 
(JTJ) with taking into account Q brings 



isTr*) ~ \w v r, v ) = -KTf^ . (4) 



An important point regarding the averaged equations (jlj is that in this form they are just 
algebraic relations between components of the smoothed hydrodynamic energy-momentum 
tensor and the average products of the metric tensor by the Ricci tensor (g ae r e p) and can- 
not therefore be taken as field equations. By splitting the products out as (g at r e p) = 
(g ae )( r ef3) + Cp where Cp is a correlation tensor, the averaged equations (jlj) become 



(9 a£ )(re(i) - \w){r»») = -KT; (hydro) - CI + \b^C\. (5) 



Here (g a/3 ) and (r aj g) denote the averaged inverse metric and the Ricci tensors which are 
supposed to describe the gravitational field due to the matter distribution T^( hydro ) _ A 
simple important observation JTJ, now is that the averaged Einstein equations (0) 
are still not "real" field equations - just a definition of the correlation tensor Cp as a 
difference between (0} and The origin of this fundamental fact is that the average 

1 For discussion on the other physical settings on general relativity facing the same problem see |14|. 

2 The mixed form is preferable here for the reason that it contains only products of metric by curvature. 
On contrary, the covariant or contravariant forms of the Einstein equations have triple products of metric 
by metric by curvature. 
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of the non-linear operator of (0) on the metric tensor g pa is not equal in general 3 to an 
operator of the same form on the average metric (g pa )'- 

( (a ae rep - ^T»>) [9 P .]) * (i9 ae )(rep) ~ ^(iTX^)) [(g P ,)] . (6) 

In order to return them the status of the field equations one must define the object 
and find its properties using information outside the Einstein equations. 

To resolve the averaging problem, and to consider it in a broader context as the 
problem of macroscopic description of gravitation, the approach of macroscopic gravity 
has been proposed 0, [TUJ , [HI, [12], [13 , C3, OH, [23 (see [II], [13] for discussion of 
the problem and references therein, J3] for discussion of the physical status of general 
relativity as either a microscopic or macroscopic theory of gravity). A covariant space- 
time volume averaging procedure for tensor fields [H], [E], [13, nas been defined and 
proved to exist on arbitrary Riemannian space-times with well-defined properties of the 
averages. Upon utilizing the averaging scheme, the macroscopic gravity approach has 
shown that (z) averaging out Cartan's structure equations brings about the structure 
equations for the averaged (macroscopic) non-Riemannian geometry and the definition 
and the properties of the correlation tensor Cp, (ii) the averaged Einstein's equations (J3 
become then the macroscopic field equations and they must be supplemented by a set of 
differential equations for the correlation tensor, (in) it is always possible to extract the 
field operator of the form for the Riemannian macroscopic metric tensor G pu and its 
Ricci tensor M pv with all other non-Riemannian correlation terms going to the right-hand 
side of (0) to give geometric correction to the averaged energy-momentum tensor 7^( hydro ) _ 
It is been also shown [TT], J2], [IS] that only in case of neglecting all correlations of the 
gravitational field the averaged equations (J3) becomes the macroscopic Einstein equations 
with a continuous matter distribution 

G°"M eP - \^G^M pv = - K T; (hydro) , (7) 

which reveals the physical status of using the standard procedure in cosmology of claiming 
(J3) to be the Einstein equations (J7J) after substitution of the matter model Q- The 
physical meaning, dynamical role and magnitude of the gravitational correlations must 
be elucidated in various physical settings. There is some evidence that they cannot be 
negligible for cosmological evolution (see, for example [23 f° r an estimation of the age of 
Universe in a second order perturbation approach). 



3 Macroscopic media in general relativity 

Derivation of the macroscopic (averaged) Maxwell field operator in macroscopic elec- 
trodynamics is easily accomplished due its linear field structure and the main problem 

3 It should noted that the inequality (jBJ has been observed in all possible averaging settings, for 
example, for a volume space-time averaging in [5], [TJ|, in the framework of a kinetic approach in [51] 
and for a statistical ensemble averaging in [22]. Relations between different averaging procedures are 
discussed in ^J, [T2] . 
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consists in the construction of models of macroscopic electromagnetic media (for exam- 
ple, diamagnetics, magnetics, waveguides, etc.) jUj, |2ZI, [28] - which relates to the 
structure of the averaged current. In general relativity the problem of construction of 
macroscopic gravitating medium models is hardly elaborated due to the following rea- 
sons: (a) existing mathematical and physical difficulties in establishing the form of the 
averaged (macroscopic) operator in (JSJ for the field equations of macroscopic gravity re- 
cedes the interest in development of macroscopic gravitating media; (b) posing on its own 
more or less realistic problem with a discrete matter creates mathematical and physical 
problems due to nonlinearity and non-trivial geometry of gravitation, to mention, for ex- 
ample, the iV-body problem, the problem of statistical description of gravity, etc. and (c) 
being relied on physically motivated phenomenological arguments (uniformity, isotropy, 
staticity, etc.) most applications of general relativity deal with effective continuous media 
if even a starting physical model is discrete in its nature like in cosmology (see Section 
EJ) or in description of extended bodies in general relativity (see jHj for discussion of the 
physical status of general relativity). 

The kinetic approach in physics is known to provide a general scheme for introduction 
of characteristics of continuous media with a known distribution function of a discrete 
configuration. But the advantages of such generality are often greatly weakened in partic- 
ular applications by difficulties of solving the Boltzmann equation to find a distribution 
function of interest. This applies to a great extent to general relativity where despite the 
formulation of the general relativistic Boltzmann equation [2H] , jHDJ the kinetic approach 
still remains useful for general definitions and considerations rather than being a working 
tool (see, for example, [21]) to derive a specific model of a macroscopic medium. 

In case of the macroscopic electrodynamics together with the volume space-time aver- 
aging on Minkowski space-time the formalism of statistical distribution functions has been 
utilized (see [27j and references therein) and it is of importance for the mathematically 
well-posed derivation of the macroscopic theory and the general structure of averaged 
current starting from microscopic electrodynamics of point-like moving charges. Further 
application of the macroscopic theory requires usually mainly phenomenological consid- 
erations to establish material relations between macroscopic average fields and induction 
field necessary to make an overdetermined system of macroscopic equations determined. A 
correct derivation of material relations is known to require [31] averaging the microscopic 
equations with a given microscopic matter model during accomplishing an averaging of 
the microscopic field equations. Though it is the only self-consistent way, the elaboration 
of such kind of approach still remains a challenge even for simple physical settings. 

On the other hand, volume (space, time, space-time) averaging procedures maintain 
their importance, direct physical meaning in application to macroscopic settings and their 
extreme clearness and descriptiveness. A volume averaging is also known to be unavoid- 
able in all macroscopic settings (including statistical approaches) [2B] 5 [HI]- EH] (HHJ and 
space-time averages of physical fields are known to have the physical meaning as directly 
measurable quantities pEI] . [35] (for discussion see [TT], [T2] and references therein). That 
greatly motivates and supports interest in applying approaches with various averaging 
schemes in physics despite corresponding (mostly mathematical, not physical) difficulties 
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in the rigorous formulation of averaging procedures. 

The paper aims to approach the problem of construction of gravitating macroscopic 
media in general relativity by using an appropriate space-time averaging scheme. 



4 Szekeres' gravitational polarization tensor 

The approach of Szekeres [3] to construct a tensor of gravitational quadruple polarization 
is based on the following assumptions: (a) the linearized theory of gravity on Minkowski 
background; (b) the linearized field equations are taken as the linearized Bianchi iden- 
tity to employ an analogy between gravitation and electromagnetism; (c) the covariant 
method of molecular moments of Kaufman is applied to construct a tensor of quadruple 
gravitational polarization. 

The equations under consideration are the contracted Bianchi identities 

^ pupa nJ 'pup (8) 

where Cp Upa is the Weyl tensor interpreted as free gravitational field, J pvp is a kind of 
"matter current" for the energy-momentum tensor t pv 

J pup — J[pu]p — ~(tp[p;u] ~ 7^9plpt,u]), (9) 

Jpu P ' p = 0. (10) 
equations (JBJ) are analogous the Maxwell equations 

fpu' U = ~Jp (11) 
C 

with (|1U|) being comparable with the conservation of the electromagnetic current j p 

V M = 0. (12) 

Let us consider a number of particles labelled by % and having masses rrii and which 
are moving in there own effective gravitational field along world lines (r,). A physical 
parameter which characterizes such a distribution is a typical characteristic distance I 
between neighboring particles. Then the corresponding microscopic energy-momentum 
tensor has the form 

m^^5 4 [x - z^r^dr,. (13) 

Assume now that due to gravitation the particles form into groups, a kind of gravitational 
molecules, which will be labelled by index a. From the physical point of view that means 
the presence of another parameter which is a characteristic size (diameter) L of such a 
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molecule with L ^> I. It is a longwave macroscopic parameter and its presence in a 
microscopic system will be defining the dynamics of the system on the distances of order 
of L. The microscopic energy-momentum tensor (|13p becomes now 



ma r c\t r\7 P r\7 v 
a i 

where y%{T a ) is a world line of the a-th molecule center of mass (3] and sf (r a ) = (rj) — 
Va( T a) is a vector connecting the i-th particle with the center of mass of the molecule 
including this particle. Let us apply now the method of molecular moments of Kaufman [3] 
to represent (fTij) as a series expansion in powers of under assumption that the effective 
gravitational field which is created by moving gravitational molecules is a weak field 
and the perturbations of the gravitational field due to relative oscillations of gravitating 
particles in molecules are small compared with the mean effective field. After averaging 
out over a typical size of the gravitational molecule (for an averaging procedure see 
[15]) one gets 4 in accordance with the Szekeres procedure 

^(molcc)^ _ rp(hec)fiu _|_ jypup ^ _|_ Q\ivpo 

where T ( - iree ^ u is the energy-momentum tensor of molecules, which has the form similar 
to (fT3j) with substitution i by a, D pup is the tensor of gravitational dipole polarization 
that can be incorporated into the quadrupole term, and Q pvpa the tensor of gravitational 
quadrupole polarization 

qp» p « = c -i { J2 I q^8\x - y a )dr a } } (16) 

which has the symmetries of the Riemann tensor. The expression for the covariant grav- 
itational quadrupole moment q a iupa is defined as 

c pa = 9»:<<-€<<-c<<+9 p :<<+ 



where 



i x ' 

O . A~- / A J 1 Ao v rl Q p rla a \ 

K vp ° = \Y. m ^r [^tt 3 * ~ f^x 1 • ( 2 °) 

3 ^ an \ dr a dr a dn dr a J 



4 No explicit averaging procedure had been used in PJ and averaged relations and equations were being 
written rather on the basis of heuristic considerations than a rigorous analysis. 
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Upon averaging (JHJ) over the typical size of gravitational molecule the following equations 
were obtained: 

(CWr = k(J^) (21) 



where 



or 



and 



= -(tfS^U + l^^U ( 22 ) 

(23) 



/ T(micro)\ _ T(frce) _ p ,cr /nt\ 
\ u pvp I ~ pup 1 livpo ■ 

The expression (I24j) is analogous to the expression for the averaged electromagnetic current 
^•(micro)^ f or a ] 3unc ] 1 Q f charged particles moving along their world lines in the effective 
electromagnetic field in accordance with the microscopic equation (jllj) with j( micro )M when 
particles are grouped into molecules jlj 

^■(micro)^ = ^(free> _ tfW ^ (25) 

where the polarization tensor P^ u is defined as an average of the quadruple polarization 
moment of the molecules p£ u (see |3] , [3] for details) 

^ = (E / dr a jC5\x-y a )). (26) 
Then equations (|21j) can be rewritten as the macroscopic equations 5 

Epupcr' — ^Jpup (27) 

for the gravitational induction tensor E^ vpa defined as 

The macroscopic equations (|27|) are analogous to the Maxwell macroscopic equations 
obtainable by means of averaging the microscopic equations (jllj) with j( micro )M with taking 
into account (J25J) 

Atc 

H ^ u = _j{^e)p ( 2 g) 

C 

for the electromagnetic induction tensor defined as 

= {D + A<kP 11v . (30) 

Unfortunately, at this point the analogy between the electromagnetism and gravitation 
which holds on the level of (jHJ), (fTU|) and (jllj) . (fT2"|) breaks. Indeed, the formal similarity 



'Gravitational macroscopic equations similar to (j21jl are known to have been proposed first in |36|. 
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of (|27j) , (J2HI an d > does not possess the structural analogy between averaged elec- 
tromagnetism and gravitation: (A) the gravitational induction tensor E^ vpa does not have 
any more the symmetries of the Weyl tensor compared with keeping the symmetries 
of f^ u ; (B) it is constructed from the second derivatives of the polarization tensor Q pvpa 
compared with the linear algebraic structure of the electromagnetic induction tensor H pv 
in terms of the polarization tensor P^ v - it is thus impossible to proceed with the formu- 
lation of phenomeno logical material relations between E pvpa and (C pupa ) as possible in 
electromagnetism (relations between H pu and (f^), or amongst the fields E, D, B, H 
and J). 

Even more important issue is that analysis of the macroscopic field equation 

(pepper) P = n(J pup "*) = Hjjivp ^ — Ppvpu j (31) 

O(e) O(l) 0(e 2 ) 

where e is a parameter measuring the value of deviation from the flat space, requires 
one to put into agreement the orders of magnitude of all quantities and reveals that the 
linearized Weyl tensor should be zero under averaging {Cp Vpo ) = 0. The Weyl tensor 
must be estimated in the perturbation theory up to the second order as it was done for 
the matter current jjj^p in the left-hand side of (jHJ). So considering some physical 
features of the polarization tensor neither the macroscopic equations (|21|). nor any other 
field equations had in fact been employed j^j. 
On the basis of the expression 

(f (micro) } = T (frcc) + l _Q ^ (32 ) 

the following material relations have been suggested 

QiOjo — (Gij)N = e g Ciojo (33) 

where N is the average number of molecules per unit volume, Gy is the quadrupole 
moment of a molecule 

Gij = / p(x)8xiSxjd 3 x, (34) 



p = p(x) is the matter density in molecules, 5xi is a vector between neighboring particles 
of a molecule. The quantity e g has been called the gravitational dielectric constant and 
in Newtonian approximation found to be 

1 mA 2 c 2 , , 

4 Uq 

where A is the average linear dimension of a typical molecule, m is the average mass of 
the molecules, Uq is a typical frequency of harmonically oscillating particles in molecules. 
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5 The Macroscopic Gravity Equations 



A macroscopic gravity approach in the perturbation theory up to the second order on 
a particular background space-time taken to be a smooth weak gravitational field is ap- 
plied to write down a system of macroscopic field equations: Isaacson's equations with a 
source incorporating the quadruple gravitational polarization tensor, Isaacson's energy- 
momentum tensor of gravitational waves and energy-momentum tensor of gravitational 
molecules and corresponding equations of motion. 

The gravitational field created by a number of particles represented by a microscopic 
energy-momentum tensor £°( micro ) j s defined by Einstein's equation 

9 at rep ~ \Vr, v = -< micro) (36) 

where k = 8nG/c 4 is Einstein's gravitational constant and G is Newton's gravitational 
constant. The Einstein equations ()36|) for the microscopic distribution of gravitational 
molecules (fTlj) have the following form: 

^-^,V%, = -< m ° lec) - (37) 

Averaging the left-hand side of the Einstein equations (}3Tj) following the Isaacson's 
high-frequency approximation approach [7J, [SJ, using the averaging procedure (one 
can also use Isaacson's averaging procedure [7J, [HJ, see also [2H]) and with taking into 
account the expression [Hj for the tensor of gravitational quadrupole polarization Q^ vpa 
in terms of the covariant gravitational quadrupole moment q^ pa (llfij) brings the averaged 
Einstein equations in the form: 

*g! - \^r {0) = -<nr ] + + l^Q^n, (38) 

where T^ W ^ is Isaacson's energy-momentum tensor of gravitational waves [7J, jH] and 
T^u^ is the energy-momentum tensor of molecules 

T$r\x) = c ~ l Y.j m ^^ 6 ^ x - ya^)]dr a . (39) 

All members in equation (J38)) can be shown to be of the same order of magnitude 0(1/ L 2 ). 
The macroscopic equations (}3*Hj) give the equations of motion for molecules 6 

T (free)^ = ^ (40) 

conservation of the energy-momentum of gravitational waves 

T (GW)^ = 0> ( 41 ) 

6 Under assumption that the background metric in the left-hand side of (|38(l represents a weak gravi- 
tational field on the flat background one can use the covariant derivatives with respect the metric in all 
relations instead of partial derivatives with respect to the flat metric. 
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and an identity for the gravitational polarization 



Q^ aix = 0. (42) 

The system of equations (|HH j) - (j4^|) is underdetermined because there are 20 unknown 
components of the tensor of gravitational polarization. It is possible to formulate two 
natural material relations. The first relation connects the traceless part of the quadrupole 
polarization tensor 

Q fipua Q ppua ^ ( QpvPpo ~\~ QpoPpv QpoPpv ~\~ QppPpa) ~\~ 

T^Sig^gpc, - g^Qpv), (43) 



where P pa = Q^ PflCT , S = P? , with the traceless energy- momentum tensor of gravitational 
waves 

c 2 ~ 

~^Qppva P — ^T^ u , (44) 

where A = X(x) the gravitational radiation polarization factor. Relation ()44|) can be 
shown to be always valid in the geometrical optics limit. 

The second material relation connects the remaining part of the polarization tensor 
Qixpvai namely its trace P pa , with a projection of the curvature tensor on the world line 
of an observer (the electric part of the curvature tensor) 

P pa = eRf pva u»u\ (45) 

where is the observer 4- velocity (4- velocity of the molecule center of mass) and e = e(x) 
is the macroscopic medium polarization factor. The relation (j45|) can be shown to lead 
to the correct expression for the 3-tensor of the average quadrupole gravitational moment 
|3] so that 

Ppu = (Poo = 0, P , = 0, Pa = (Gij)N). (46) 

where N is the average number of molecules per unit volume, (Gij) is the averaged 
quadrupole moment of a molecule (}3*4*j) . Then the material relation 3j can be recovered 
in the form 

QiOjo — (Gij)N = e g R i0 j (47) 

that gives e = e g with the gravitational dielectric constant e g defined in jH] as (|35|). 

Thus the system of equations (|HH |) - (j4T|) . (j33j), (j43|) is fully determined and can be used 
to find the gravitational and polarization fields for the macroscopic gravitating systems. 



6 The Static Weak-field Macroscopic Medium 

The averaged Einstein equations (|3*K|) have been derived under assumption of the weak 
gravitational field though they can be considered to be formally valid for any background 
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metric g$ with a given tensor of gravitational quadrupole polarization Q^ upa and the 
material relations f)44j) and (|4"H1 ). The definition of the tensor of gravitational quadrupole 
polarization Q^ upa (llfjjl adopted here is essentially valid on the flat space-time background 
due to the used definitions of molecule's center of mass [3], [HZ]. Therefore, the averaged 
Einstein equations (|3*Hjl with ()16|) can be only consistently applied in the framework of 
the weak-field approximation. 

A model of static weak-field macroscopic medium with quadrupole gravitational po- 
larization is considered here. The model is based on three assumptions. 

(1) Newtonian gravity conditions for the energy-momentum tensor of molecules T$v, 

rp({ree) ^(free) rp (free) ,7-, (free) ,7-, (free) _ rpifrea) u, v _ 2 ( aq\ 

1 00 1 ij 1 1 00 1 Oi 1 1 00 — 1 iiv U U — , 

where the observer 4-velocity (4-velocity of the molecule center of mass) u u is u v = 
(1, 0, 0, 0) and \x is the macroscopic matter density. The condition (]4*%|) means that the 
gravitational field created by gravitational molecules is essentially Newtonian. 

(2) The macroscopic metric tensor gf} is static, 

m = °' (49) 

which means that is the macroscopic matter density fi depends only on spatial coordinates, 
ji = (J>{x a ). 

(3) The condition of the weak field approximation for the macroscopic metric tensor 

(o) 

9$ = V»v + e V' ( 50 ) 

where is rj^ the flat space-time metric, 7/00 = — 1, r}n = 1, ^22 = 1, ^33 = 1 and rj^ = 
if \i 7^ v, h^ u is the arbitrary perturbation functions depending here only on spatial 
coordinates and e is the smallness parameter, e<C 1. 

In general relativity the conditions (|48|) -(|5U| ) for a microscopic energy- momentum ten- 
sor 7^( micro ) anc i f or the metric tensor g 00 = — (l + fr) , 9oi — 0, g%j = 0, are known to lead 
to the Newtonian limit of the Einstein equations (J36j) to result in the Poisson equation 
for the Newtonian gravitational potential if = (p(x a ) 

A(p = AnGfx. (51) 

The Newtonian limit of the macroscopic gravity equations (|3*H1 ) under conditions 
-([50)1 for the macroscopic tensor g^J 



^ = -(l + f), <?[? = !, 9^ = 1, $=1, 9® = 0, (52) 

should bring a generalization of the Poisson equation which incorporates the effect of 
gravitational quadrupole polarization. 
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For the case of static weak-field macroscopic medium Isaacson's energy-momentum 
tensor of gravitational waves Tji„ w ^ vanishes 

?^ W) = 0, (53) 

and no gravitational radiation polarization factor A is involved 

A = 0. (54) 

The the macroscopic medium polarization factor e is here 

(55) 

with the gravitational dielectric constant e g defined |3] as (|33j) . 



7 The Equation for the Macroscopic Gravitational 
Potential 

Calculation of the equation for the macroscopic gravitational potential ip from the macro- 
scopic gravity equations (J3%|) under conditions (JH^-flEDI for the macroscopic tensor g$ 
([52]) brings the equation 

Aip = AttG/2 + AV (56) 

3<r 

where A 2 y? = A(A^) is the Laplacian of the Laplacian of ip. This is a non-trivial gener- 
alization of the Poisson equation for the gravitational potential (p of Newtonian gravity 
(|5ip. This is a biharmonic equation due the presence of the term A 2 (p. The equation (|56|) 
involves a singular perturbation, since in case of the vanishing gravitational dielectric con- 
stant, e g = 0, this equation becomes the Poisson equation, but if e g ^ 0, this equations 
change its operator structure to be of the fourth order equation in partial derivatives of 
ip as compared with the Poisson second order partial differential equation. 
It is convenient to introduce the factor 

1 4vrGe n 



k 2 " 3c 2 



(57) 



with k having a physical dimension of inverse length, [k 2 ] = length 2 . Then the equation 
f[56[) takes the form 

Ay? = 4ttG7i + AY (58) 

By using the definitions of the gravitational dielectric constant e g (J33J), the charac- 
teristic oscillation frequency of molecule's constituents Uq, macroscopic matter density 
fi = 3m/4nA 3 and the average number of molecules per unit volume N = 4irD 3 /3 with D 
as a mean distance between molecules, the factor k~ 2 can be shown to have the following 
form 

7E\-Fs) A - ( 59 ) 



k 2 46 \D 3 J 
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Here the dimensionless factor 9, 

* - A* m 

reflects the nature of field responsible for bounding of discrete matter constituents into 
molecules. If 9 ~ 1, the molecules of self-gravitating macroscopic medium are considered 
to be gravitationally bound. For instance, considering a macroscopic model of galaxy 
as a self-gravitating macroscopic medium consisting of gravitational molecules taken as 
double stars, 9 m 1 as such galactic molecules are gravitationally bound. If one takes the 
molecules to be of electron-proton type, like atoms, the factor 9 « 10 40 , which makes the 
factor k~ 2 essentially insignificant. 

The dimensionless ratio A/D reflects the structure of macroscopic medium. If (A/D) « 
1, the macroscopic medium behaves itself like a liquid or solid. If (A/D) < 1, the macro- 
scopic medium behaves itself like a gas. For the macroscopic galactic model for the present 
epoch the macroscopic medium is like a gas, since (A/D) ~ 10 _1 — 10~ 2 , which makes 
the factor A 3 /D 3 to be of order of 10~ 3 — 10~ 6 . However, for earlier times of galaxy 
formulation this factor can be expected to be of much greater order of magnitude up to 
1 - 10. 

It is useful to introduce a dimensionless factor k~ 2 L~ 2 , 

1 1 /A 3 \ fA 2 \ , , 

^ ™ T? . (61) 



k 2 L 2 49 \D 3 J \L 2 J 

where L is the characteristic scale of the macroscopic gravitational field. The dimension- 
less ratio A/L reflects the scale of significant change in the macroscopic gravitational field. 
If (A/L) ~ 1, the macroscopic gravitational field changes significantly on the scale L. If 
(A/D) <C 1, the scale L does not reflect the presence of gravitational molecules. 

Thus, the structure of the factor k~ 2 in (|59p. or k~ 2 L~ 2 in (|6ip. is model depen- 
dent, and its particular value is fully determined by a particular model of self-gravitating 
macroscopic matter. When a macroscopic medium has the factor 

1 1 /4 3 \ /A 2 ' 

" S ' 1, (62) 



k 2 L 2 49 \D 3 J \L 2 

the equation for the macroscopic gravitational potential <p is (158(1 which is convenient to 
write in the following form 

AV - k 2 Ap = -4nGk 2 fi. (63) 
It can be rewritten as a system of two second order partial differential equations 

A<p = f, (64) 

Af - k 2 f = -47rGk 2 fi, (65) 
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for the unknowns tp(x,y,z) and f(x,y,z) with given fi(x,y,z) and k 2 . In this case of 
the significant gravitational quadrupole polarization of a macroscopic medium, which will 
be referred to as Case I, the equations (|63j) is singular with respect to the gravitational 
dielectric constant e g because k 2 ~ l/e 9 and the limit t g — > cannot be accomplished in 
a solution to (|54*)l and (JoT)|) . Equations (JMJ) an d (jSH)) do not have a limit to the Poisson 
equation (J51|) of Newtonian gravity, nor does a solution to (|64|) and (|65|) have a limit to a 
solution to (joTj) . It should be pointed out that the equations (EU) and (J63j) are essentially 
valid either inside the macroscopic matter source, or outside the matter source in the 
close vicinity of the macroscopic matter configuration boundary where the effect of strong 
gravitational quadrupole polarization is still significant, 

~ ^Av? » Ay?. (66) 

As far as asymptotically Aip/k 2 L 2 becomes much less than Aip the equations ()64)1 and 
(jHSJ) are not valid anymore. See the general solution to equations (|6~4"|) and (|63|) and their 
asymptotic behavior in Section |H| 

When a macroscopic medium has the factor 

^hi Ti «L (67) 



k 2 L 2 46\D 3 J \L 2 

the equation (JH%|) for the macroscopic gravitational potential i/j can be rewritten effectively 
in the form of the Poisson equation 

= 4ttG> (l + J^) , (68) 

for the unknowns if>(x, y, z) with given n(x, y, z) and k 2 L 2 . It is possible because the A 2 ijj 
term in (|H%1 ) is well approximated as 

^A 2 ij w ^A^ < A^. (69) 

In this case of the insignificant gravitational quadrupole polarization of a macroscopic 
medium, which will be referred to as Case II, the equation (jQHl) is nonsingular with 
respect to the gravitational dielectric constant e g because 1/ k 2 ~ e 5 and the limit e g — > 
can be accomplished in the solution to ()68p. Equation does have a limit to the 
Poisson equation (JSTj) of Newtonian gravity, and a solution to (p)8*j) always has a limit to 
a solution to (f5lf . It should be pointed out that the equation (J68|) holds either inside 
the macroscopic matter source, or outside the matter source. A solution to equation (jHEJ) 
gives the macroscopic gravitational potential ijj(x,y,z) defined everywhere. 

The equation (|55|) describes also asymptotic behavior, L 2 — > oo, for the field of 
the macroscopic gravitational potential ip(x, y, z)for a macroscopic matter distribution 
with strong gravitational quadrupole polarization. A solution to equations (JMj) and (|65|) 
must be matched to a solution to equation (|55|) on a surface outside the source between 
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the source boundary and the nearby zone to get the macroscopic gravitational potential 
ip(x,y,z) in asymptotic region and, as a result, to define it everywhere. 

The Cases I and II have different equations and different physics, which is reflected 
by solutions to equations (|64p and and equation 



8 Case I: The Strong Quadrupole Polarization 

The system of equations (JMj) and (|65J) can be solved by separation of variables. Consider 
the equation (|fi5j) in the spherical coordinates (x, y, z) — > (r, 8, <fi) when the Laplacian 
becomes 

&--—(r 2 —\- — L 2 -— coth— 1 92 (7Q) 
r 2 dr \ dr J r 2 ' d8 2 88 sin 2 8 d(f) 2 

Applying the method of separation of variables for the function f(r, 8, <p) in with 
assuming the angular part is given by the spherical harmonics Y^d, <fi) [3*%j . 

/(r,M) = — (71) 
r 



one can represent the Green function in a spherical harmonics expansion [39| as 

00 1 

G(n,r 2 ) = J2J2 giiru^Y^^Yl^fo), (72) 

;=0 m=-l 

which gives the equation for the radial Green function gi(r%, r 2 ) 
d 2 S 



r i-^j[ r i9i( r i,r 2 )] - k 2 r 2 gi(r l ,r 2 ) - 1(1 + l)gi(r l ,r 2 ) = -An5(ri - r 2 ). (73) 
It has the solution 

gi{ri,r 2 ) = kiiikr^k^kr^, ij(Ar), r x < r 2 , k t (kr), n > r 2 (74) 
and the solution for f(r, 8, <fi) for r\ > r 2 has the form 



00 1 

/( ri ) = Gk 2 ^ Yl ^i{^i)Y l m {e U (j)i) x 

1=0 m=-l 

J / u(r 2 )^(A;r 2 )F r | l (^,0 2 )rf0 2 sin^^ 2 r 2 rfr 2 . (75) 

This is a multipole expansion of where the particular structure of the function of 

a macroscopic matter distribution /i(r 2 ) will make the particular structure of multipole 
expansion. 
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Now one can solve the Poisson equation ()64j) with known /(r) to obtain ip(r). But 
rather than continuing solving the system in spherical harmonic representation, let us 
solve the system in general form by using Green's function method. The solutions are 
very useful for illustration of the character of the macroscopic gravitational potential 
<p(r) for the static weak- field approximation with the quadrupole moment tensor. One 
can always accomplish expansion of the general solution with respect to the spherical 
harmonics to get a multipole expansion. 

The system of equations (|54*|) and (jo^j) must be supplemented by boundary conditions 
for the unknown /(r) and (p(r) with given /i(r) and k 2 : 

lim <f(r) = exists and is bounded, lim <p(r) = 0, < |r| < oo, (76) 

| r| — ^0+ |r|^oo 



lim /(r) = exists and is bounded, lim /(r) = 0, < |r| < oo. (77) 

|r|^0+ |r|^oo 

The system of partial differential equations (JMj) and (|65|) with the boundary conditions 
(|76|) and (|77|) is a Dirichlet boundary value problem on interval < |r| < oo. One can 
consider also two Dirichlet problems for intervals (0, |r |) and (|r | , oo) where r is the 
radius vector of a macroscopic matter configuration. 

To solve first the nonhomogeneous modified Helmholtz equation (165(1 . one needs to 
find the corresponding Green function 

AG/(n, r 2 ) - fc 2 G/(n, r 2 ) = -AttS^ - r 2 ) (78) 

with the boundary condition 

lim G?,(n,r a )=0. (79) 

|ri|-»oo 

The Green function can be found 1 39 1 to be 



„-fc|ri-r2 

G / (r 1 ,r 2 ) = - (80) 

|ri - r 2 | 

and the solution for /(r) is 

r v(r 2 )dV 2 . (81) 

|ri - r 2 | 

The second equation (j51|) has the Green function 

G y ( ri ,r 2 )= 1 (82) 
|ri - r 2 | 

as the general solution to the Green equation, 

AG^(r 1 ,r 2 ) = -4vr5(r 1 -r 2 ), (83) 
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to bring the solution for </?(r) 

An J |ri-r 2 | 

Now the solution for the equation (jUHj) can be written as 

^ = "IT J J |r 2 -r 3 | ^ (85) 

If the macroscopic gravitational potential <p(r) is calculated in the origin of the coordinate 
system, that is ri = 0, r 2 = R, the formula (fKKj) becomes 

Gk 2 fdV R f dV r e- k ^- r \ 
V = -^IW\I |R~r| Mr) ' (86) 

The asymptotic form of the solution (jSfij) as |R| — > + is given by 

lim ld ~ -4vrGA; 2 u(r) |R| 4 . (87) 

|r|-»0+ 

It shows that the macroscopic gravitational potential (p (r) has better analytical properties 
as |r| — > + then the gravitational potential of the Poisson equation. 
The asymptotic form of the solution (jHUjl as |R| — > oo is given by 



(1 R R p- fe l R l 

lim p ~ -Gk 2 -^——^ (88) 

|r[-»«> |R| 3 

where Qij is the total quadrupole moment of the macroscopic mass distribution 



Qij = J fi(r)r irj dV r . (89) 

The formula ()88|) shows that the macroscopic gravitational potential (p(r) fades out in the 
vicinity of the macroscopic matter source boundary on the characteristic distance k~ x . 
For the distances |R| ^> k^ 1 the inequality does not hold and the inequality (JoT?|) 
is valid. Therefore to find a proper asymptotic form of the macroscopic gravitational 
potential f(r) for |R| ^> k^ 1 one should solve the equation (|ES|) and match the solution 
(jnni)) see Sectional with the solution (|HU|) . 

9 Case II: The Weak Quadrupole Polarization 

The equation (J68j) can be solved by the Green function method. To solve the Poisson 
equation (p5)l. one needs to find the corresponding Green function 

AG^( ri ,r 2 ) = -47r<f(ri-r 2 ) (90) 
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with the boundary condition 



lim G^(n,r 2 )=0. (91) 

|ri|— too 



The Green function can be found |39| to be 



G^ ri ,r 2 ) = - r (92) 

I"l — f2 



and the solution for ip(r) is 



*M = - a (^m?)I\A K - (93) 

If the macroscopic gravitational potential ip(r) is calculated in the origin of the coordinate 
system, that is ri = 0, r 2 = R, the formula (J85j) becomes 

Behavior near the center and at infinity are essentially the same up to renormalization 
by the factor (1 + l/k 2 L 2 ). The multipole expansion jlU] for the macroscopic gravitational 
potential ip with L ~ Rq has the form 

^- G (^ik){^wM + -}- (95) 

The factor (1 + l/k 2 L 2 ) renormalizes the whole potential to contribute into all terms in 
the multipole expansion. 
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